Dynamical gap generation in graphene nanoribbons: An effective relativistic field 

theoretical model 
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We show that the assumption of a nontrivial zero band gap for a graphene sheet within an 
effective relativistic field theoretical model description of interacting Dirac electrons on the surface 
of graphene describes the experimental band gap of graphene nanoribbons for a wide range of widths. 
The graphene band gap is dynamically generated, corresponding to a nontrivial gapless solution, 
found in the limit of an infinitely wide graphene ribbon. The nanoribbon band gap is determined 
by the experimental graphene work function. 

PACS numbers: 72.80.Vp,11.10.Kk 
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I. INTRODUCTION 

Graphene sheets, carbon nanotubes (CNTs) and 
graphene nanoribbons (GNRs) are very closely related 
subjects governed by quantum mechanical effects in a 
constrained tvsro-dimensional (2D) surface. In this com- 
pactified dimensional world, symmetry and scale drive 
the astonishing properties that these one-atom-thick ma- 
terials exhibit. GNRs are recognized today as promising 
nanosheets to be used in nanoelectronic and spintronic 
devices^. Such material has attracted a lot of atten- 
tion since free-standing graphene sheets were reported 
by using different experimental techniques. Experiments 
confirmed that the charge carriers can be described as 
massless Dirac fermions^, opening a solid basis for its rel- 
ativistic treatment. It is expected that massless fermions 
moving through graphene have an approximate ballis- 
tic transport behavior with very small resistance due 
backscattering suppression^. Graphene is also a good 
thermal conductor, when compared with other semicon- 
ductor junctions, due to the high mobility of carriers at 
room temperaturo^i^. 

Energy band gaps for GNRs weremeasured in Refsi^ 
audi using different methods. First principles ah initio 
methods^, density functional theory (DFT)^, and tight 
binding (TB) models^^ were applied to calculate GNR 
band gaps for a reasonable range of not too small widths. 
In opposition to TB models, first-principles calculations 
do not support metallic nanoribbons^. On the other 
hand DFT calculations predict energy gap oscillations 
as a function of the GNR for very small widths^. 

The two-dimensional (2D) graphene sheet is essentially 
a zerogap semiconductor (with the Fermi level Ep pre- 
cisely at = 0), with a linear chiral carrier kinetic energy 
dispersion relation given hy E = fivp fc^, where k is the 
2D carrier wave vector and Ui? is the Fermi velocity (in- 
dependent of carrier density) . 

The fermion behavior on the graphene sheet has been 
discussed within the free hamiltonian. Ho, picture using 



four different approaches^i: (i) as Schrodinger fermions 
Ho = /2m* , where m* is an effective mass; (ii) as 
ultra-relativistic Dirac fermions, where Ho = cap; (iii) as 
massless Dirac fermions with Ho = vpap, where a are 
the Pauli matrices; and (iv) as massive chiral fermions 
with Ho = a p^ jlm* , where ct is a pseudospin matrix de- 
scribing the two sublattices of the honeycomb lattices^. 

This work is based on a relativistic effective field the- 
ory of electrons containing a local four-fermion interac- 
tion. The effective interaction Lagrangian parameterizes 
the actual forces acting on the electrons. It resumes the 
relevant dynamics from the atom-electron and electron- 
electron interaction within the many-body system (in- 
cluding e.g. electron-phonon interactions). The ratio- 
nale to use effective contact interactions goes back to the 
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that the Hartree-Fock approximation of the 



single particle energy of the many-fermion system with 
the point-coupling interaction can be compared to the 
first order contribution of the surface density in a den- 
sity functional theoryi^. 

The relativistic picture of interacting fermions has to 
be understood together with the concept of an effective 
mass. In the case of the graphene the effective mass 
vanishes, corresponding to a gapless material. A non- 
vanishing effective mass for GNRs means that a gap is 
open. Our model provides a mechanism for dynamical 
mass generation for the nanoribbons, even for a gapless 
graphene. 

We assume an homogeneous 2D system of electrons (or 
holes),— where the free part of the electron Hamiltonian 
corresponds to case (iii), i.e.. Ho ~ vpSp. The two sub- 
lattices of graphene, associated to the Dirac points K and 
K', are taken into account by introducing a factor 2 in 
the electron multiplicity. 

The model is based on a relativistic Lagrangian field 
theoretical approach in which a four-fermion scalar and 
a vector point-like interactions are considered and the 
many-fermion problem solved in the Hartree-Fock ap- 
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FIG. 1: Schematic graphene bands for conducting (upper) 
and valence (lower) fermions. For the discussion of A see the 
text. 



being the largest integer smaller than A VF/tt. The max- 
imum momentum of the negative energy electrons filling 
the valence band, as illustrated in Fig. [U is A. 

The fermion mass, i.e. the electron self-energy, is gen- 
erated dynamically and is computed from the gap equa- 
tions for the GNRs, using a self-consistent Schwinger- 
Dyson equation at one-loop level (Hartree-Fock approxi- 
mation). Recall that the Schwinger-Dyson equation has 
a cutoff in the momentum of the negative electron en- 
ergy in the valence band given by A. The gap equations 
for the self-energies, after integration on the longitudinal 
momentum, arc given by 



proximation on a 2D surface. The model has only two 
free parameters, one related to the strengths of the four- 
fermion interaction and the Fermi momentum cut-off, A. 
Exploring the connection of the model with graphene, 
i.e., by assuming that graphene is a nontrivial gapless 
solution of the dynamical equations, we are able to re- 
late the strengths of the interactions and cut-off with the 
graphene work function. This procedure allows to fix the 
parameters and calculate the gaps of the nanoribbons 
without any further assumptions. 

Our work is organized as follows. In section II, we 
present the model and derive the gap equations. In sec- 
tion III, we present the solution of the gap equations for 
the nanoribbon and compare our results with the exper- 
imental datai^ Our concluding remarks are given in the 
last section. 



where 



^(E 
2GA? 



2G 



S{k) = In 



A2-fc2 



E 



A2 - fc2 



(2) 

(3) 
(4) 

(5) 
(6) 



II. THE MODEL 

The electron-hole dynamics is described by a relativis- 
tic quantum field effective Lagrangian in a 2-t-l space- 
time, where the interaction has a pairwise and contact 
form, which generalizes a nonrelativistic 2D modcU^, 

C = §(17^9^ - m)* + Gs(**)2 + G„(§7^*)2, (1) 

where s and v refer to scalar and vector couplings, respec- 
tively, and refers to fermion fields. In the following, we 
use a system of units such that h ^ vp = 1. 

Our method of solving the model for the nanoribbon 
starts from a nanotube geometry where the fermions are 
constrained to move on the cylindrical surface. The ap- 
propriate coordinates are the position along the nanotube 
symmetry axis (z) and the angle in the transverse plane 
(0) - see Ref^ for details. The GNRs are obtained by 
unfolding an infinite length carbon nanotube. The re- 
sulting sheet is a GNR with a width W, after eliminat- 
ing dangling bonds. The procedure to obtain GNRs from 
CNTs is described in Refi^. 

The Dirac modes are quantized along the transversal 
direction after demanding that the spinors vanishes at the 
edges of the nanoribbon. From this condition it follows 
that the allowed fermion transverse momenta are given 
by kn = mr/W, where n = ±1, ±2, • 



and a2 (fc„ - E«)2 + (S, + m)^. 

In Eqs. da]) to dll) we define x = AW/tt, G = - 3 
and a factor of 2 was included to take into account for 
electrons at the K and K' points. Note that the scalar Gs 
and the vector Gy coupling constants combine into a sin- 
gle coupling constant G. Furthermore, the self-energies 
given by Eqs. (Il])-(|3]) include a rest fermion mass m. 
However, hereafter we consider only the m = case. 

The scalar self-energy is identified with the half 
bandgap nannoribbon. The self-energy vector part is 
actually a quadri-vector containing space-time compo- 
nents. However, due to symmetry only the time compo- 
nent given by E° survives. This component is directly 
related to the number of particles N through 
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(7) 



where L is the longitudinal length of the GNR. 

In a rclativistic picture, the infinite graphene sheet 
{W — >■ oo) has no electrons in the conducting band. The 
valence band is fully occupied and there is no gap be- 
tween both bands, i.e. Ss = 0. Then, it remains to 
parametrize only Eq. ([3]), since Ss = is a trivial solu- 
tion for Eq. ([2]). Recall that we are considering massless 
fermions m = 0. Anyway, insofar considering very large 
but finite values for W, ([5]) admits other solutions besides 
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the trivial = and one has an opening gap generated 
dynamically. 

Let us look for the nontrivial solution of ([2]). This 
defines a unique coupling constant, a function of A, G = 
Gcrit, which will be called the critical graphene limit. If 
one takes the continuum limit of Eq. ([2]) with m ~ 0, 
after setting y = fc„/A it follows that 
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Given that the above integral is equal to tt/2, then the 
critical graphene limit requires a 
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The A dependence can be eliminated in favor of the 
graphene experimental work function value WF= 4.8 
eV— . The work function is related with the chemical 
potential fi, which needs to be computed in a consistent 
thermodynamic The pressure P, the density 

energy £ of the system determine the chemical potential 
by = {P + £)/cr, with cr = N/A and where A is 
the area of the GNR. P and £ are given by the diagonal 
terms of the energy-momentum tensor constructed from 
the model Lagrangian (jlj in the usual way. In a mean 
field approacbi^, 
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where n(x) = VA^ — x, and 



P = < T^^ > 
_ -1 



(11) 



where z is the longitudinal direction of the GNR. Then, 
the expression for the chemical potential reads 



GA2 



fix) - VA^ + Sf = -WF, 



(12) 



where WF is the work function. As a side remark, ac- 
cording to referencei^, the chemical potential computed 
in the way described above, contains already both the 
bulk and surface contributions from the energy density 
and pressure. 

Graphene is recovered in the limit — > oo and, in this 
case, X = AW/n also diverges, while f{x — > oo) = tt/2. 
Therefore, for graphene if one sets G = Gcrit and = 0, 
one defines the parameters of the model. Then, using 
Gcrit, it follows from Eq. ^ that A = 2WF. Given that 
the experimental value for the graphene work function 



is 4.8 eV— with determine that A = 9.6 eV. From 
now on, we will use always these results for computing 
the nanoribbon gaps within the effective relativistic field 
theoretical model given by the lagrangian density ([T} . 

The above definitions fully parameterizes Eqs. (Hl-Q, 
allowing the computation of a self-consistent solution for 
Eg as a function of the nanoribbon width W. Further- 
more, given a value for A, one can compute the surface 
density cr = A^/tt . 



III. RESULTS 

In order to solve the model numerically, besides Gcrit 
and A = 9.6 eV, we use vp = c/300fiS, Our numerical 
procedure is as follows: i) For a given W, x = AW/tt 
and k„ = wk/W are defined, constraining the integer 
n-max < X] a) given that E^ = 0, fc^ + E^ and Eq. 

^ can be solved iteratively to compute Eg = 2Es. E° 
is calculated directly from Eq. (j3]). 

In Fig. [2] the results of the model for the band gap 
Eg as a function of the nanoribbon width are compared 
with the experimental data of reference^. As observed 
in the figure, the model reproduces the GNRs experi- 
mental gap for a large range of nanoribbon widths. We 
call the reader attention, that the curve in Fig. [2] is the 
direct outcome from the model without any ad hoc nor- 
malization to reproduce the right Eg scale. This strik- 
ing agreement between experimental data and theoretical 
predictions seems to indicate that the model, as defined 
above, captures the essential physics required to under- 
stand the gap formation in nanoribbons. Moreover, the 
results summarized in Fig. [2] give support to the idea 
that effective interacting valence Dirac electrons which 
generates dynamically a band gap is useful to describe 
the electronic properties of GNRs. 

The experiments reported in Refi^ andii show that 
for CNTs with radius between 0.5 and 1.5 nm, the work 
function oscillates with an amplitude of ^ 0.5 eV or less 
around the graphene work function WF. As the CNT 
radius increases, the oscillations rapidly become very 
smooth. This consideration may help in understanding 
why our model still works for GNRs with small widths. 
Moreover, it also shows that the graphene scale, used in 
the definition of our model, is robust regarding the GNR 
finite-size details. Possibly, the massless fermions picture 
reinforces this point, since the localization of carriers is 
strongly suppressed favoring ballistic behavior for them. 

Let us now compare the model with other theoreti- 
cal predictions. In Refj^, first-principles calculations for 
three types of armchair-edged GNRs have been reported. 
Their results are similar to ours. For smaller widths, pre- 
vious theoretical calculations^ have been done using DFT 
functionals from the works of Perdew et al. (PBE)^ and 
Heyd et al. (IISE)^^. According to the authors, the 
bandgap oscillates as a function of W . A similar behav- 
ior is observed in our model - see Fig. |31 Typically, our 
results are in between those of the two DFT (PBE, HSE) 
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the physics require to understand the GNR from smaU 
to large widths. Moreover, the graphene scale is used to 
defined the model parameters and it seems to be a ro- 
bust scale regarding the GNR finite-size details. In this 
sense, our results are very auspicious for further investi- 
gations of GNRs with the present model which may, as 
claimed before, represent a first order approximation to 
a DFT calculation. A nice feature of the relativistic field 
approach is its simplicity. Indeed, it allows for a detailed 
analytical description that, eventually, can guide our in- 
tuition in the investigation of GNRs and in looking for 
general driving physics trends. 




FIG. 3: GNR bandgap as a function of ribbon width. Theo- 
retical calculations from PBE and HSE density functionals^. 



calculations. 

From the point of view of the relativistic model, the os- 
cillations in the physical quantities are related to rimax- 
As the width W decreases, x = AW/tt decreases and 
higher values for the integer n become possible. When- 
ever Umax increases by unity, due to the behavior of the 
derivative of N relative to the single particle energy, a 
Van-Hove singularity shows up. This explains the oscil- 
latory behavior of the GNR bandgaps observed in our 
model. As seen in Fig. |31 whenever W gets larger, the 
amplitude of oscillations diminishes and they become siz- 
able only for a width of the order of a few nanometers or 
smaller. 



IV. CONCLUDING REMARKS 

Our results for the effective relativistic field theoretical 
model show that it accounts reasonably well for some of 



To conclude, we would like to make general considera- 
tions on how to improve the present effective model. For 
small widths, due to the increase on the kinetic energy, 
some of the fermions can jump into the conduction band. 
Therefore, for small W's the model instead of having a 
single cut-off A, will require two, one for the valence band 
Ai and another for the conduction band A2. The pres- 
ence of the new scale could lead to larger amplitude oscil- 
lations as the ab initio calculations suggests - see figure 
[31 Moreover, when W becomes smaller, the edges distin- 
guish armchair and zigzag GNRs, as discussed in Refi^. 
This requires refined boundary conditions for the spinors 
at the GNR at the boundaries. In general vanishes 
due to the symmetry between fc„ and — fc„. However, 
if this symmetry is broken due to the boundary condi- 
tions, then one should solve Eqs.([2])-(l4]) self-consistently. 
In any case, the relativistic model still keeps its simplic- 
ity and the way we have solved the model is, essentially, 
unchanged. 
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